Effective band-structure in the insulating phase versus 
strong dynamical correlations in metallic VO2 
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Using a general analytical continuation scheme for cluster dynamical mean field calculations, 
we analyze real-frequency self-energies, momentum-resolved spectral functions, and one-particle 
excitations of the metallic and insulating phases of VO2. While for the former dynamical correlations 
and lifetime effects prevent a description in terms of quasi-particles, the excitations of the latter 
allow for an effective band-structure. We construct an orbital-dependent, but static one-particle 
potential that reproduces the full many-body spectrum. Yet, the ground state is well beyond a 
static one-particle description. The emerging picture gives a non-trivial answer to the decade-old 
question of the nature of the insulator, which we characterize as a "many-body Peierls" state. 

PACS numbers: 71.27. +a, 71.30.+h, 71.15.Ap 



Describing electronic correlations is a challenge for 
modern condensed matter physics. While weak corre- 
lations slightly modify quasi-particle states, by broad- 
ening them with lifetime effects and shifting their ener- 
gies, strong enough correlations can entirely invalidate 
the band picture by inducing a Mott insulating state. 

In a half-filled one-band model, an insulator is re- 
alized above a critical ratio of interaction to band- 
width. Though more complex scenarios exist in realis- 
tic multi-band cases, a common feature of compounds 
that undergo a metal-insulator transition (MIT) upon 
the change of an external parameter, such as temper- 
ature or pressure, is that the respective insulator feels 
stronger correlations than the metal, since it is precisely 
their enhancement that drives the system insulating. 

In this paper we discuss a material where this rule of 
thumb is inverted : We argue that in VO2 it is the insu- 
lator that is less correlated, in the sense that band-like 
excitations are better defined and have longer lifetimes 
than in the metal. Albeit, neither phase is well described 
by standard band-structure techniques. Using an an- 
alytical continuation scheme for quantum Monte Carlo 
solutions to Dynamical Mean Field Theory (DMFT) Q, 
we discuss quasi-particle lifetimes, k-resolvcd spectra (for 
comparison with future angle resolved photocmission ex- 
periments) and effective band-structures. While dynam- 
ical effects are crucial in the metal, the excitations of 
the insulator are well described within a static picture : 
For the insulator we devise an effective one-particle po- 
tential that captures the interacting excitation spectrum. 
Still, the corresponding ground state is far from a Slater 
determinant, leading us to introduce the concept of a 
"many-body Peierls" insulator. 

The MIT of VO2 has intrigued solid state physicists 
for decades @, i, i, H, @, 0, & S, EE EE El EE 13 ■ A 
high temperature metallic rutile (R) phase transforms at 
T c =340 K into an insulating monoclinic structure (Ml), 
in which vanadium atoms pair up to form tilted dimers 



along the c-axis. The resistivity jumps up by two orders 
of magnitude, yet no local moments form. Despite exten- 
sive efforts, the mechanism of the transition is still under 

debate @,0,S,SE3,Eil,E3- two scenarios compete : In 



the Peierls picture the structural aspect (unit-cell dou- 
bling) causes the MIT, while in the Mott picture local 
correlations predominate. 

VO2 has a d 1 configuration and the crystal field splits 
the 3d-manifold into ti g and empty e a g components. The 
former further split into e" g and a\ g orbitals, which 
overlap in R-VO2, accounting for the metallic charac- 
ter. Still, the quasi-particle peak seen in photocmission 

(pes) @, Ed, Eii 

is much narrower than the Kohn-Sham 
spectrum of density functional theory (DFT) in the local 
density approximation (LDA) 0], and eminent satellite 
features evidenced in PES are absent. In MI-VO2, the 
a\ g form bonding/antibonding orbitals, due to the dimcr- 
ization. As discussed by Goodenough [3j , this also pushes 
up the e g relative to the a\ g . Yet, the LDA Q yields a 
metal. Non-local correlations beyond LDA were shown 
to be essential EE E3- Indeed, recent Cluster DMFT 
(CDMFT) calculations [la], in which a two-site vanadium 
dimer constituted the DMFT impurity, opened a gap, 
agreeing well with PES and x-ray experiments ll|, [l2| . 

Starting from these LDA + CDMFT results (la] for the 
Matsubara t2 g Green's function G(tuj n ) we deduce the 
real frequency Green's function G{ui) by the maximum 
entropy method [HI and a Kramers-Kronig transform. 
The self-energy matrix E(w) we obtain by numerical in- 
version of G(tu) = J^k [ w + I 1 ~ -^k — ^(^P 1 Hi w itb 
the LDA Hamiltonian H, and the chemical potential /x. 

Fig. [1] shows (a) the diagonal elements of the R- 
VO2 self-energy, and (b) the resulting k-resolved spec- 
trum. Notwithstanding minor details, the a\ g and 
self-energies exhibit a similar dynamical behavior. The 
real-parts at zero energy, 0?S(O), entailing relative shifts 
of quasi-particle bands, are almost equal, congruent with 
the low changes in their occupations vis-a-vis LDA [ID], 
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FIG. 1: (color online) Rutile VO2 : (a) self-energy (S — fx). 
Real (imaginary) parts are solid (dashed). As comparison Ml 
SyE ai , Sffi 01 g -aig are shown, (b) spectral function A(k,w) 
and solutions of the QPE (blue). The LHB is the (yellow) 
region at -1.7 eV, the broad UHB appears (yellow) at ~2.5 eV. 



and with the isotropy evidenced in experiment [19( . 

Neglecting lifetime effects (i.e. QEwO) , one-particle ex- 
citations are given by the poles of G(oS) : det[wk + /-* — 
i/k — 5ffi(tJk)] = 0- We shall refer to this as the quasi- 
particle equation (QPE) [23)]. For static or absent this 
reduces to a simple eigenvalue problem. In regions of low 
the QPE solutions will give an accurate description 
of the position of spectral weight and constitute an effec- 
tive band-structure of the interacting system. Yet, due 
to the frequency dependence, the number of solutions is 
no longer bounded to the number of orbitals. 

Below (above) -0.5 (0.2) eV, the imaginary parts of 
the self-energy - the inverse lifetime - of R-VO2 is con- 
siderable. Due to our limited precision for 3E(0), we 
have not attempted a temperature dependent study to 
assess the experimental bad metal behavior, but the re- 
sistivity exceeding the Ioffc-Regel-Mott limit [2C| indi- 
cates that even close to the Fermi level, coherence is not 
fully reached. At low energy, the QPE solutions (dots 
in Fig. [TJd) closely follow the spectral weight. Above 
0.2 eV, regions of high intensity appear, howbeit, the 
larger broadens the excitations, and no coherent fea- 
tures emerge, though the positions of some e g derived 
excitations are discernible. At high energies, positive and 
negative, distinctive features appear in 3£(cj) that are 
responsible for lower (upper) Hubbard bands (L/UHB), 
seen in the spectrum at around -1.7 (2.5) eV. The UHB 
exhibits a pole-structure that reminds of the low-energy 
quasi-particle band-structure. Hence, an effective band 
picture is limited to the close vicinity of the Fermi level, 
and R-VO2 has to be considered as a strongly correlated 
metal (the weight of the quasi-particle peak is of the or- 
der of 0.6). This is experimentally corroborated by the 
fact that an increase in the lattice spacing by Nb-doping 
results in a Mott insulator of rutile structure Q ■ 

The imaginary parts of the Ml a\ g on-site, and a\ g - 
a\ g intra-dimer self-energies, Fig. [Ij,, are larger than 



in R-VO2, usually a hallmark of increased correla- 
tions. However, we shall argue that correlations are in 
fact weaker than in the metal. Indeed, the dimeriza- 
tion in Ml leads to strong inter-site fluctuations, evi- 
denced by the significant intra-dimer a\ g ~a\ g self-energy. 
Fig. [2] displays the MI-VO2 self-energy in the a\ g bond- 



ing/antibonding (bab) basis, E 



b/ab 



The a\ g (anti)bonding imaginary part is low and varies 





FIG. 2: (color online) Self-energy (E-fj.) of MI-VO2 in the a lg 
bab-basis : (a) real parts. The black stripes delimit the a lg 
LDA bandwidths, dashed horizontal lines indicate the values 
of the static potential A. (b) imaginary parts. Self-energy 
elements are dotted in regions irrelevant for the spectrum. 



little with frequency in the (un)occupied part of the 
spectrum, thus allowing for coherent weight. In the 
opposite regions, the imaginary parts reach huge val- 
ues. The e g elements are flat, and their imaginary parts 
tiny. This is a direct consequence of the drastically re- 
duced e g occupancy which drops to merely 0.14. These 
almost empty orbitals feel only weak correlations, and 
sharp bands are expected at all energies. A first idea 
for the a\ g excitations is obtained from the intersections 
u)+fx— e(,/ O 6(k)=0?Sb/ a (,(w) as depicted in Fig. [2^, where 
the black stripes delimit the LDA a\ g bandwidths. The 
(anti)bonding band appears as the crossing of the (blue) 
red solid line with the stripe at (positive) negative en- 
ergy. Hence, the (anti)bonding band emerges at (2.5) 
—0.75 eV. Still, the antibonding band is much broadened 
since 3S a 6 reaches -1 cV. To confirm this, we solved the 
QPE and calculated the k-resolved spectrum (Fig. |3Ja,). 
As expected, reasonably coherent weight appears over 
nearly the entire spectrum from -1 to +2 cV, whose po- 
sition coincides with the QPE poles : The filled bands 
correspond to the a\ g bonding orbitals, while above the 
gap, the e g bands give rise to sharp features. The anti- 
bonding ai g is not clearly distinguished since e g weight 
prevails in this range. The L/UHB have faded : a mere 
shoulder at -1.5 eV reminds of the LHB. Finally, con- 
trary to R-VO2, the number of poles equals the orbital 
dimension. Since, moreover, the real-parts of the Ml- 
VO2 self-energy are almost constant for relevant ener- 
gies 24J , we construct a static potential, A, by evaluat- 
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ing the dynamical self-energy at the LDA band centers 
(pole energies) for the e g (ai g ), sec Fig. [2^ 0. Fig. [3)3 
shows the band-structure of H^+A : The agreement with 
the DMFT poles is excellent. Our one-particle potential, 
albeit static, depends on the orbital, and is thus non- 
local. We emphasize the conceptual difference to the 
Kohn-Sham (KS) potential of DFT : The latter gener- 
ates an effective one-particle problem with the ground 
state density of the true system. The KS energies and 
states are auxiliary quantities. Our one-particle poten- 
tial, A, on the contrary, was designed to reproduce the 
interacting excitations. The eigenvalues of H^+A are 
thus not artificial. Still, like in DFT, the eigenstates are 
SDs by construction, although the true states are not (see 
below). The crucial point for MI-VO2 is that spectral 
properties are capturable with this effective one-particle 
description. It is in this sense that MI-VO2 exhibits only 
weak correlation effects. The weight of the bonding ex- 
citation is Z=(l — d w $lY,ij(uj))~=_ 7aV ~0.75, and thus 
larger than the rutile quasi-particle weight (see above). 

What is at the origin of this overall surprising coher- 
ence? For the e g orbitals, this simply owes to their deple- 
tion. For the nearly half-filled a\ g orbitals the situation is 
more intricate. It is a joint effect of charge transfer into 
the a\ g bands, and the bonding/antibonding-splitting. 
Indeed, the filled bonding band experiences only weak 
fluctuations, due to its separation of several eV from the 
antibonding one. To substantiate these qualitative argu- 
ments, we resort to the following model, which treats the 
solid as a collection of Hubbard dimers : 

H = -*£ { C lla C l2* + h - C )~ t ^2 C ia C lHa + U "«t n KJ. 
la «=1,2 il 

(1,1') 

(1) 

Here, c\ ia (cu a ) creates (destroys) an electron with 
spin a on site i of the I th dimer. t is the intra-dimcr, 
t± the inter-dimer hopping, U the on-site Coulomb 
repulsion, and we assume half-filling. First, we discuss 
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FIG. 3: (color online) MI-VO2 : (a) spectral function A(k,^). 
(blue) dots ((a) & (b)) are solutions of the QPE. (b) The (red) 
dots are the eigenvalues of i/k+A. See text for discussion. 



the t±^ limit, which is an isolated dimer : the 
Hubbard molecule. We choose t—0.7 eV, the LDA 
intra-dimcr ai g -ai g hopping, and U =4.0 eV [HI for 
all evaluations. The bonding/antibonding-splitting, 
A bab =-2t + \/16i 2 + U 2 =3A8 eV, gets enhanced with 
respect to the U=0 case. In MI-VO2, the embedding 
into the solid, and the hybridization with the e g reduce 
the splitting to ~3 eV, as can be inferred from the one- 
particle poles (Fig. [3]), consistent with experiment [ill ]. 
The ground state of the dimer is given by \ipo) = 
{At/(c-U)(\l T> — I T l» + (ITI 0) + |0 nm/aM 
which is intermediate to the Slater determinant (SD) 
(the four states having equal weight), and the Heitler- 
London (HL) limit (double occupancies projected out). 
With the VO2 parameters, the model dimer is close 
to the HL limit @. The inset of Fig. gb shows the 
projections of the ground state onto the SD and the 
HL state. The former, | (SD\ipo) \ 2 , equals the weight 
of the band-derived features in the spectrum (for U>0 
satellites appear), while the other measures the double 

occupancy J2i( n i1 n il) = 1 — K-^-^IV'o)! • For [7=4.0 eV 
the latter is largely suppressed, as a consequence of the 
interaction : The N-particle state is clearly not a SD. 
Still, the overlap with the SD, and thus the coherent 
weight, remains significant, i.e. one-particle excitations 
survive and lifetimes are large. To do justice to the 
seemingly opposing tendencies of correlation driven 
non-SD-bchavior, coexisting with a band-like spectrum, 
we introduce the notion of a "many-body Peierls" state. 

The charge transfer from the e g into the then almost 
half-filled a\ g orbitals, finds its origin in the effective re- 
duction of the local interaction in the bab-configuration : 
While for U=A eV, (SD\H\SD) = 2.0 eV in the SD limit, 
it reduces to merely (ipo\H\tpo) = 0.91 eV in the ground 
state. In fact, inter-site fluctuations are an efficient way 
to avoid the on-site Coulomb repulsion. In MI-VO2, this 
effect manifests itself in a close cancellation of the local 
and inter-site self-energies in the (un-)occupied parts of 
the spectrum for the (anti)bonding a\ g orbitals. 

The gap-opening in VO2 thus owes to two effects : The 
sclf-cncrgy enhancement of the a\ g bab-splitting, and 
a charge transfer from the e g orbitals. The difference 
in corresponds to this depopulation, seen in exper- 
iments and theoretical studies d, , and leads to 
the separation of the a\ g and e g at the Fermi level. The 
local interactions thus amplify Goodenough's scenario. 

To show that the embedding of the dimer into the 
solid docs not qualitatively alter our picture of the Ml 
phase, we solve the model, Eq. (fT]), using CDMFT. This 
moreover allows to study the essentials of the rutile to 
Ml MIT by scanning through the degree of dimeriza- 
tion t at constant interaction strength U and embedding, 
or inter-dimer hopping, tj_. For the latter we assume 
a semi-circular density of states D±(lu) of bandwidth 
W=4:t±. In MI-VO2, the t± for direct a\ g -a\ g hopping 
is rather small, yet e^-hybridizations lead to an effective 
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Dj^-bandwidth of about 1 eV. We choose U=4t±, and an 
inverse temperature (3=10/t±. Fig. [4^ displays the or- 
bital traced local spectral function A(uj)=A b (uj)+A ab (uj) 
(b,ab denoting again the bonding/antibonding combi- 
nations) and the bonding self-energy X b (w) for differ- 
ent intra-dimer hoppings t : In the absence of t, the 
result equals by construction the single site DMFT so- 
lution (£ b =E ab ), which, for our parameters, is a corre- 
lated metal, analog to R-VO2. The spectral weight at the 
Fermi level is given by A b/ab (0) = D ± (±t - 3ffi b/ab (0)), 
with ^.T, b/ab (0)=T^S ab {0). Thus a MIT occurs at 
t + 5RS ab (0) = 2t±, when all spectral weight has been 
shifted out of the bandwidth : Above t/t±=0.5 we find a 
many-body Peierls phase corresponding to MI-VO2. In 
Fig. 2k we have indicated again the graphical QPE ap- 
proach : The system evolves from three solutions per or- 
bital (Kondo resonance, L/UHB) at t=0 to a single one at 
t/t±=Q.6. Hence the peaks in the insulator are not Hub- 
bard satellites, but just shifted bands. The embedding, 
t± , broadens the excitations and washes out the satellites 
of the isolated dimer, like for MI-VO2. Still, as a function 
of t, the coherence of the spectrum increases, since the 
imaginary part of the (anti-)bonding self-energy subsides 
at the renormalized (anti-)bonding excitation energies. 
Our model thus captures the essence of the rutile to Ml 
transition, reproducing both, the dimerization induced 
increase in coherence, and the shifting of excitations. 
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Peierls singlets as in Ml, while the tilted pairs are in a 
Mott state. Hence, we interpret the seminal work of [1] 
as the observation of a Mott to many-body Peierls tran- 
sition taking place on the tilted chains when going from 
M2 to Ml. To illustrate this, we solve again Eq. |T]) for 
appropriate parameters. The tilted M2 chains are akin 
to the rutile phase, yet with a reduced a\ g bandwidth 0]. 
Thus we now choose U=6t±, f3=10/t±, and vary t. All 
solutions shown in Fig. [4Jd are insulating, however, the 
diverging self-energy at vanishing intra-dimer coupling 
(t=0, tilted "M2" chains) becomes regularized with the 
bond enhancement (t>0, "Ml"). The imaginary part of 
the self-energy gets flatter and the system thus more co- 
herent. The above is consistent with the finding of (S=0) 
S=l/2 for the (dimerized) tilted pairs in M2-V0 2 

While our results do not exclude surprises in the direct 
vicinity of T c 



22] 



the nature of insulating VO2 is shown 
to be rather "band-like" in the above sense. Our analyti- 
cal continuation scheme allowed us to explicitly calculate 
this band-structure. The latter can also be derived from 
a static one-particle potential. Yet, this does not im- 
ply a one-particle picture for quantities other than the 
spectrum. Above all, the ground state is not a Slater de- 
terminant. Hence, we qualify MI-VO2 as a "many-body 
Peierls" phase. We argue that the weakness of lifetime 
effects results from strong inter-site fluctuations that cir- 
cumvent local interactions in an otherwise strongly cor- 
related solid. This is in striking contrast to the strong 
dynamical correlations in the metal, which is dominated 
by important lifetime effects and incoherent features. 
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point. We thank AIST, Tsukuba, for hospitality. JMT 
was supported by a JSPS fellowship. Computer time was 
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FIG. 4: (color online) (a) spectral function (top), real 
(middle), imaginary (bottom) bonding self-energy Ej, of the 
CDMFT solution to Eq. (TJ for f/=4.0ix, /3=10/t ± , and 
varying intra-dimer hopping t/t±. 5RS(,(o;)=— 5RE a t(— u), 
9E; ) (cij)=Q ! E £l i,(— uj) by symmetry, (b) Imaginary Matsuba- 
ra self-energy, Q'Eb(iui)=QY, a b(iu) , for U=6t±, /3=10/t± and 
varying t. Inset: Projection of the SD and HL limit on the 
Hubbard molecule ground state (t—0.7 eV, t±=0) versus U. 



Under uni-axial pressure or Cr-doping, V0 2 develops 
the insulating M2 phase Q in which every second vana- 
dium chain along the c-axis consists of untiltcd dimers, 
whereas in the others only the tilting occurs. We may 
now speculate that the dimerized pairs in M2 form a\ g 
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